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ABSTRACT 

Context. Large grids of synthetic spectra covering a widespread range of stellar parameters are mandatory for different stellar and 
(extra-)Galactic physics applications. For instance, such large grids can be used for the automatic parametrisation of stellar spectra 
such as that performed within the AMBRE project for which the main goal is the stellar atmospheric parameters determination for the 
few hundreds of thousands of archived spectra of four ESO spectrographs. 

Aims. To fulfil the needs of AMBRE and future similar projects, we have computed a grid of synthetic spectra over the whole optical 
domain for cool to very cool stars of any luminosity (from dwarfs to supergiants) with metallicities varying from 10~ 5 to 10 times the 
solar metallicity, and considering large variations in the chemical content of the cr-elements. 

Methods. For these spectrum computations, new generation MARCS model atmospheres and the Turbospectrum code for radiative 
transfer have been used. We have also taken into account as complete as possible atomic and molecular linelists and adopted, in the 
spectral synthesis, the same physical assumptions and input data as in the MARCS models. This allowed us to present a grid with a 
high consistency between the atmosphere models and the synthetic spectra. 

Results. A new grid of 16 783 high resolution spectra over the wavelength range 3 000 to 12000 A has been computed with a 
spectral resolution always larger than 150000. Normalised and absolute flux versions are available over a wide range of stellar 
atmospheric parameters for stars of FGKM spectral types. The covered parameters are 2 500 K<r eff <8 000 K, -0.5<log(g)<5.5 dex 
and -5.0<[M/H]<+1.0 dex. Moreover, for each combination of these stellar parameters, five different values of the enrichment in 
cc-elements have been considered (0.0, ±0.2 dex and ±0.4 dex around the standard values). This library is thus relevant for any stellar 
type and luminosity class, present in old and intermediate-age stellar populations from extremely metal-poor to metal-rich in content. 
This grid is made publicly available through the POLLUX database (about 50% of the spectra are already included in this database) 
and in FITS format upon request to the authors. 

Key words, stars: atmospheres - stars: abundances - stars: fundamental parameters - surveys - astronomical databases: miscellaneous 



1. Introduction 

AMBRE is a project defined between the Observatoire de la 
Cote d'Azur and the European Southern Observatory (ESO). Its 
main goal is the automated stellar parametrisation of the FEROS, 
HARPS, UVES and GIRAFFE archived spectra. This parametri- 
sation of more than 300 000 spectra will be made publicly avail- 
able through the ESO archive. Furthermore, AMBRE offers the 
opportunity to test on a substantial amount of real spectra the 
performances of parametrisation algorithms such as MATISSE 
dRecio-Blanco et al. 20061 or DEGAS ( |Kordopatis et al. 201 l| l 
which have been specifically developed for the analysis of 
the Radial Velocity Spectrometer (RVS) data collected by the 
ESA/Gaia mission. We refer the reader to |Worley et al. 2012 for 
a detailed description of AMBRE and the complete analysis of 
the FEROS spectra. 

Within the AMBRE project, the adopted algorithms for the 
automatic parametrisation of stellar spectra rely on a library of 
reference spectra to derive the atmospheric parameters and the 
chemical abundances. Preparing such a spectra grid, covering a 
wide range of parameters is mandatory in order to allow the rapid 
analysis with adapted algorithms of large quantities of observed 
spectra from any stellar type, such as those provided by current 



(RAVE, SEGUE, Gaia/ESO Survey.. ) and future (Gaia/RVS) 
spectroscopic surveys. Therefore, when planning AMBRE, li- 
braries of observed (empirical) and/or computed (synthetic) stel- 
lar spectra were investigated. We ultimately concluded that there 
is the lack of an extensive library of stellar spectra with the fol- 
lowing properties: (i) coverage of a wide range of atmospheric 
parameters, (ii) coverage of a very large spectral domain in the 
optical, (iii) high spectral resolution, (iv) well-established atmo- 
spheric parameters, (v) inclusion of chemical abundance varia- 
tions. 

Indeed, although important progress has been made over the 
last few years to build high-quality empirical libraries of stellar 
spectra (see, in particular, the last release of the CFLIB library by 
IWu et al. 201 lb . observed spectra libraries still often suffer from 
severe shortcomings. For instance, most of them rely on rather 
low resolution spectra (CFLIB or MILES. IVazdekis et al. 20101 1 
or rather restricted wavelength coverage (as the ELODIE li- 
brary, Prugniel & Soubi ran 200 l| l and/or suffer from flux cali- 
bration uncertainties. Furthermore, the atmospheric parameters 
attached to these compiled stellar spectra may be subject to cau- 
tion due to the lack of determination of some parameters, differ- 
ent adopted bibliographic sources and/or method of determina- 
tion (see Soub iran et al. 2 010). However we point out the impor- 
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tant effort conducted recently by I Wu et al. 20 1 1 1 to obtain a ho- 
mogeneous set of stellar atmospheric parameters. But this new 
CFLIB library still suffers from a rather incomplete coverage 
of the parameter space (particularly for cool and/or metal-poor 
stars). Furthermore, no information on the chemical abundance 
of a elements is provided (to our knowledge, all available em- 
pirical libraries provide only the main three, in the best cases, at- 
mospheric parameters, i.e. the effective temperature, the surface 
gravity and the mean metallicity). All of these reasons prevented 
the use of these empirical libraries within AMBRE because, as 
mentioned previously, the main goal of AMBRE is to derive, 
from high-resolution spectra covering very different regions of 
the whole optical domain, the atmospheric parameters together 
with the [a/Fe] enrichment of the target stars for the purposes of 
galactic archaeology. 

Regarding the synthetic spectra libraries that are avail- 
able, the situation was also not completely satisfactory. 
Mun ari et al. 20051 present an extensive grid of cool to very 
hot stars based on Kurucz/ATLAS9 model atmospheres and 
the SYNTHE spectral synthesis code (see also their Tab. 1 for 
a comparison of the principal characteristics of most publicly 
available libraries before 2005). Although the number of avail- 
able spectra is impressive (because the authors considered dif- 
ferent values of microturbulent and rotational velocities, spec- 
tral resolution and sampling law), these spectra actually corre- 
spond to just a few thousand combinations of T e g , log(g), [M/H] 
values. Furthermore, the largest spectral resolution {R ~20000) 
is too low for the analysis of the ESO spectra and, as for 
the empirical libraries, no variations of the [a/Fe] enrichment 
have been considered (two values of [a/Fe] were computed at 
three distinct metallicities, only). Finally, the metallicity range is 
also too small for general parametrisation (-2.5 < [M/H]< 0.5). 
Another synthetic high-resolution spectra grid has been pre- 
sented by Coe lho et al. 20 05 again based on Kurucz/ATLAS9 
models. This grid covers the whole optical domain but, once 
again, the adopted variations in [a/Fe] are still not satisfactory 
(only [a/Fe] = and 0.4 dex are considered for every models) for 
the automatic parametrisation purposes for which this chemical 
index is searched. Finally, [Kirby 20lT] published a grid very sim- 
ilar to our needs in terms of atmospheric parameter ranges but, 
unfortunately, it is dedicated to the red part of the optical domain 
(A > 6 300 A). These spectra are computed from a new grid 
of Kurucz/ATLAS9 atmosphere models and the MOOG code. 
There is however a lack of molecular lines considered in the 
spectral synthesis. Only the CN, C2 and MgH molecules were 
considered (partly because the number of TiO lines to be con- 
sidered led to prohibitive computational times). Therefore, the 
provided spectra for cool and/or metal-rich stars might not be 
well modelled in this grid. 

For all these reasons, we decided to compute a new synthetic 
grid of stellar spectra in order to fulfil the AMBRE needs. The 
goal of this paper is to present and make such a grid publicly 
available. The computed spectra cover the whole optical spectral 
range (from 3 000 to 12 000 A) at high resolution (R > 150 000) 
in order to be used for the analysis of most spectrograph data. 
They correspond to FGKM stars of any luminosity (except 
White Dwarfs) with a very wide range of metallicities and a- 
element abundance variations at any metallicity. Furthermore, it 
is the first grid based exclusively on MARCS model atmospheres 
and associated tools to compute the synthetic spectra (with there- 
fore the same physical assumptions and ingredients, such as the 
extensive molecular linelists). So far, it is also the most extensive 
grid by about one order of magnitude in the number of available 



spectra of cool stars (mostly because we span cooler tempera- 
tures, lower metallicities and several variations in [a/Fe]). This 
grid can easily be used as reference for the automatic determi- 
nation of atmospheric parameters (as carried out in the AMBRE 
project) or for building composite synthetic stellar populations 
for extragalactic studies. In particular, the simulation of the inte- 
grated spectra of dwarf galaxies as those found in the Milky Way 
vicinity with stars having large or small [a/Fe] values (Frebel et 
al. 2010; Letarte et al. 2010) requires such a library with a broad 
range of [a/Fe] chemical enrichments. This is also mandatory 
for future and present galactic spectroscopic surveys looking for 
accreted halo stars. 

This paper is organised as follow: Sect. 2 is devoted to the 
description of the methods, assumptions and input data adopted 
for generating the stellar spectra. The AMBRE grid of synthetic 
spectra is then presented in Sect. 3. Sect. 4 & 5 discuss some 
properties and limitations of this grid. We present in Sect. 6 
how these spectra can be retrieved. The conclusions are given 
in Sect. [7] 



2. Computation of the synthetic spectra 

The AMBRE project is primarily devoted to the analysis 
of FGKM stars, i.e. the cool stars that are the main tar- 
gets of the MARCS model atmospheres for almost 40 years 
(IGustafsson et al. 19751 >. Therefore, the presented grid is based 
on the last generation of MARCS model atmospheres presented 
in IGustafsson et al. 20081 The synthetic spectra have been gen- 
erated in a consistent way (as far as possible) with these mod- 
els. The spectra were indeed computed with the Turbospectrum 
code for spectral synthesis V9.2 (Alva rez & Plez 1998| with sub- 
sequent improvements through the years by B. Plez). This code 
uses the same routines and input data as the MARCS code. 
We remind that MARCS models are computed assuming one- 
dimensional plane-parallel or spherical geometry (depending on 
the gravity, see below), hydrostatic and local thermodynamic 
equilibria and convection is treated following the mixing-length 
theory. 

We have seen that, so far, all the other available grids of 
synthetic spectra are based on Kurucz/ATLAS9 model atmo- 
spheres. Gustafsson et a l. 20081 present a comparison of their 
new MARCS models with the ATLAS9 ones (see, also, Plez 
2011). They report a rather good agreement in the tempera- 
ture structures for solar-type or cooler dwarf stars of different 
metallicities. However, larger discrepancies appear at lower T e ff 
and/or lower [M/H]. Furthermore, for giants, MARCS models 
appear slightly cooler in the upper layers (probably due to the 
improved molecular data collected by the MARCS group over 
the years). It is therefore expected that the presented AMBRE 
spectra grid may represent an important improvement towards 
the modelling of cool stars of any luminosity (from dwarfs to 
supergiants) and metallicity. 

2. 1 . Atomic and molecular line opacities 

In the spectrum calculations a number of specific line lists for 
hydrogen, metals, and molecules covering the 3000-12000 A 
spectral range were considered. 

Atomic line data were recovered from the Vienna Atomic 
Line Database (content of august 2009; |Kupka et al. 199"9] >. Only 
lines from neutral or singly ionised atoms were kept (which is 
a reasonable assumption for the considered effective tempera- 
tures). This led to a total of about 215 000 atomic lines from 
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Table 1. Isotopic compositions adopted in the calculations of the molecular equilibrium (see, for references, the Turbospectrum 
code). 
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76 neutral metals and 60 singly ionised metals (among those 
lines, about 35 100 of them are transitions of Fe I and 33 600 
of Fe II). Due to the huge amount of adopted lines over the 
wide spectral range, no correction or astrophysical calibration 
of the atomic data were performed (the completeness and the 
quality of the VALD linelists could be questionned although this 
database is continuously improved thanks to huge efforts con- 
ducted by several groups). Therefore, we warn future users of 
this synthetic spectra grid that it is not optimised for fine anal- 
ysis or abundance determination from fitting individual spectral 
lines. As in IGustafsson et al. 20081 hydrogen self broadening is 
estimated following Barklem et al. (2000a) and damping con- 
stants for many strong metallic lines are calculated using re- 
cent data (Barklem et al. (2000b) and related references J] Other 
lines are broadened with the classical Unsold theory with an as- 
sumed enhancement factor for the Van der Waals broadening set 
to 2.5 except for Nal, Sil, Cal and Fel for which the enhance- 
ment factor was set to 2.0, 1.3, 1.8 and 1.4, respectively (see 
IGustafsson et al. 2008b . 

The considered molecular lines come from the CH, NH, OH, 
MgH, SiH, CaH, FeH, C2, CN, TiO, VO, and ZrO molecules 
with their corresponding isotopic variations. In these calcula- 
tions, we disregarded CO lines since they appear beyond the 
IR wavelength limit of the present grid at A > 1.5 microns. 
We however point out that this molecule together with sev- 
eral other molecular species have been considered in the cal- 
culation of the chemical equilibrium (see the complete list in 
Gustafsso n et al. 2008b . In total, 41 different species (see Tab[T] 
for the adopted isotopic ratios) consisting of more than 20 mil- 
lions molecular lines were taken into account (half of these 
lines being TiO transitions and one fourth belong to the ZrO 
molecule). These molecular line lists have been compiled by B. 
Plez and are identical to those used by Gustafsson et al. 2008 
(see their Tab. 2 for a list of references). We point out that, for 
molecular lines, the damping was not taken into account for the 
same reasons as in the MARCS models. 

2.2. Selection of MARCS model atmospheres 

We recovered 3 358 MARCS model atmospheres (version of 
September 2009; only converged models were kept). The two 
main atmospheric parameters of the selected models (see Tab.|2| 
are the effective temperature (T e ff) and the stellar surface gravity 



1 One can expect some differences (particularly in the H-line profiles) 
in our computed fluxes when compared to previous grids since such 
improvements have never been considered up to now by other groups 
computing large spectra grids. 



(log(g)). The effective temperatures vary between 2 500 K and 
8 000 K (step of 200 K below 4 000 K and 250 K above) and the 
logarithm of stellar surface gravities from -0.5 to 5.5 dex (step 
of 0.5 dex), g in units of cm/s 2 . 

Then, the selected models cover the mean metallicity ([M/H] 
hereafter, except the a-elements) range from -5.0 to +1.0, span- 
ning 15 different values (see Fig. [T] & |2j. We point out that, 
throughout this paper and the AMBRE project, the mean metal- 
licity refers to solar scaled abundances of all elements heavier 
than He. It should also be noted that the adopted metallicity 
range allows for analysis of very different types of Galactic stel- 
lar population, including extremely metal poor halo stars. 

Finally, a selection has been performed on the enrichment in 
or-elements. As in IGustafsso n et al. 20081 the chemical species 

0, Ne, Mg, Si, S, Ar, Ca, and Ti are considered as a-elements 
and we consider that they vary in lockstep. We remind the reader 
that [Gustafsson et al. 20081 have shown that the effects of an a- 
enhancement on the temperature structure in the atmosphere be- 
come significant below r e jj< ~4500 K for dwarfs and giants 
(leading to larger temperature gradients in the deep layers and to 
higher surface heating, both being mainly caused by TiO lines). 
For the present AMBRE grid, at a given metallicity, we only 
selected the MARCS models with an [a/Fe] enrichment consis- 
tent with the chemical properties of most of the galactic stars, 

1. e. [a/Fe]= 0.0 for [M/H]> 0.0, [ar/Fe]= +0.4 for [M/H]< -1.0 
and [a/Fe]= -0.4 x [M/H]for -1.0 < [M/H]< 0.0. 

Regarding the assumed geometry in these calculations 
(model atmospheres and synthetic spectra), plane-parallel mod- 
els have been considered for +3.5 < log(g)< +5.5 (in these cases, 
the atmospheric extension is negligible relative to the radius, and 
thus sphericity effects can be neglected). These models have a 
microturbulent- velocity parameter set to 1.0 km/s. A spherical 
geometry has been considered when log(g)< 3.0 since sphericity 
effects may be important for low masses and/or low log(g) (see 
IHeiter & Eriksson~2 006 and Gustafss on et al. 200 81. For these 
gravities, a mass of 1.0 M Q and a microturbulence parameter 
of 2.0 km/s have been considered. Fig. [T] presents the distribu- 
tion of these selected model atmospheres with respect to their 
metallicity and the adopted plane-parallel/spherical geometry. 

On the other hand, we emphasize that all the selected 
MARCS models have been computed assuming a mixing length 
parameter a = l/H p — 1.5 (Gustafss on et al. 20 08) and solar 
abundances from Grevesse et al. (2007) that reflect recent 
3D hydrodynamic, LTE-departures and improved atomic and 
molecular data for the solar modelling. 
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Table 2. Atmospheric parameters of the AMBRE grid. The description for the other parameters adopted in the spectra calculation is 
detailed in Sect. 2 & 3. We also recall that (i) for each set of (r e ff, log(g) [M/H]) values, 5 spectra with different [a/Fe] enrichments 
have been computed, and that (ii) some models are missing (approach to the Eddington flux limit and/or poor convergence). 



Atmospheric Parameter (units) Minimum Maximum Step Number of possible values 



r eff (in K) 


2500 


8000 


200 
250 


for T eff < 3900 
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for [M/H]< -3.0 
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[cc/Fe] (in dex) 


-0.4 


+0.8 


0.2 
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Fig. 1. Distribution of the selected MARCS model atmospheres with respect to the metallicity for effective temperatures cooler (left 
panel) or hotter (right panel) than 4 000 K. Filled and blank histograms refer to the adopted geometry for the radiative transfer, i.e. 
plane-parallel or spherical, respectively. 



In total, 3 358 MARCS models were selected (2 115 models 
with r eff >4000 K and 1 243 models with r eff <4000 K). We point 
out that a few models are missing (all the combinations of the at- 
mospheric parameters were not available) due to the approach to 
the Eddington flux limit and/or poor convergence of the MARCS 
code (see Gustafsso n et al. 2008b . The ranges of the atmospheric 
parameters of the AMBRE grid are summarised in Tab.|2]and il- 
lustrated in Fig.|2](see also Fig.[T|. 



3. The AMBRE grid of synthetic spectra 

From the atmosphere models presented in the previous section, 
we computed the final grid that consists of 16783 flux nor- 
malised and absolute-flux spectra (10575 spectra have r e ff> 
4000 K and 6208 have T eS < 4000 K; it has to be pointed out 
that 7 cool star spectra with available model atmospheres are 
missing due to numerical issues). For reasons of consistency 
with the MARCS model atmospheres of Gustafs son et al. 2008 
the spectra have been computed assuming the solar chemical 



composition of Grevesse et al. (2007) and similar geometry, 
mass and microturbulence parameter (see previous Section). 

In the spectrum calculations, five different [a/Fe] enrich- 
ments were adopted for any selected MARCS model. The 
adopted [a/Fe] values vary from -0.4 to +0.4 dex (step of 
0.2 dex) with respect to the canonical values that correspond 
to the original abundances of the selected MARCS models (see 
Tab. |2]and Fig. [2] bottom left panel). Therefore, these computa- 
tions were performed with chemical abundances that are some- 
times inconsistent with those adopted in the calculation of the 
model atmosphere structure. However, for each synthetic spec- 
trum, all continuous and line opacities were recomputed from 
the model structure with the right chemical composition (see 
Subsect. [5] for an estimation of the errors induced by this as- 
sumption). 

The delivered spectra cover the wavelength range between 
3 000 and 12000 A with a constant wavelength step of 0.01 A 
(900000 pixels in total). With the Nyquist-Shannon sampling 
condition, the adopted spectral sampling corresponds to an 
equivalent resolving power (R = A/AA) larger than 150000 in 
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Fig. 2. Distribution of the AMBRE synthetic spectra grid in the atmospheric parameters and [a/Fe] space. We reiterate that only 
one value of [a/Fe] for every [M/H] was adopted during the model atmosphere selection process. This is illustrated in the bottom 
right panel where the AMBRE spectra computed from MARCS models with consistent [a/Fe] ratios are plotted with filled circles 
while open circles refer to all the other AMBRE spectra computed with atmosphere models having inconsistent [a/Fe] ratios. 



the ultraviolet/blue spectral domain to ~ 600 000 in the near in- 
frared region. We point out that we decided not to perform any 
convolution (spectral resolution and instrumental profile, stellar 
rotation, macroturbulence,...) for these spectra in order to keep 
an easy and fast adaptation to the properties of any spectrograph 
and/or astrophysical application. The adopted spectral domain 
and sampling indeed cover the wavelength ranges explored by 
all the ESO optical spectrographs, which are the main targets 



of the AMBRE project. This spectra grid is also useful for al- 
most any other astronomical optical spectrograph since spectra 
with any lower resolving power (or large-scale broadening) can 
be easily computed from the provided grid. However, we warn 
the future users of this grid that the adopted wavelength step 
might not be well adapted for very high-spectral resolution ap- 
plications. Indeed, the narrowest lines could be slightly under- 
sampled in the present grid. This has been tested by computing 
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a few spectra with a wavelength step ten times smaller than the 
one of the AMBRE grid (i.e. 0.001 A) and then by degrading 
these spectra to different spectral resolutions. Small differences 
between these spectra having a better sampling and the AMBRE 
ones are found for spectral resolution larger than 50 000 (they 
can reach 2-3% in relative flux in the blue when R — 100 000 or 
higher). 

All computations were performed in serial mode on the 
952 cores (2.4Ghz Opteron) of the high-performance comput- 
ing facility SIGAMM of the Observatoire de la Cote d'Azur. 
The total CPU time used by this project was about 10 5 hours 
(~11 years if one would have used a single CPU!). This large 
CPU time is almost completely caused by the computation of 
the opacities induced by the ~ 10 7 molecular lines in the differ- 
ent atmospheric layers, and it explains why these lines have not 
been considered in most of the synthetic spectra grids available 
so far. 

Finally, the synthetic spectra are given in absolute surface 
flux at each wavelength bin (in units of erg/cm 2 /s'/A) and in 
normalised flux, relative to the local continuum flux (no units). 
The continuum-normalised spectrum is derived by dividing the 
absolute flux spectrum by the computed continuum flux in ev- 
ery spectral bin. The size of one spectrum file in FITS format is 
3.5 Mbytes (about 10 times more in ASCII format) and the total 
grid of 16783 spectra has a size of 59 Gbytes in FITS. 

4. Contribution of the different chemical species to 
the spectrum flux 

In order to study the contribution of the different chemical 
species to the total opacity, we computed test spectra by 
considering separately the different linelists. This analysis 
has been performed on the spectra of different types of stars 
that are slightly enriched in metals ([M/H]= +0.5) in order 
to enhance the different contributions. We show in Fig. [3] 
the line blocking by different species (atoms and molecules). 
The atomic lines are the dominant source of line opacity for 
dwarf and giant stars hotter than ~4 000 K. Molecules become 
the strongest contributors for cooler stars. More precisely, 
carbon-based diatomic molecules (C2, CN and CH) have an 
important contribution with respect to other molecules for 
T e g> 4 500 K. However, at lower temperatures for dwarf stars, 
the molecular contribution is then dominated by molecules 
composed of or-elements (TiO, ZrO, VO, MgH, CaH, SiH 
and OH). Their blanketing effect thus strongly increases at 
such low effective temperatures. These molecules (and TiO 
in particular) may therefore strongly affect the temperature 
structure in the atmosphere. On the contrary, we notice that cool 
giant spectra are more dominated by carbon-rich molecules 
than their dwarf counterparts with similar effective temperatures. 

We have also examined the line blocking of atomic and 
molecular lines for two typical giant stars of the galactic bulge 
and halo (Fig. |4] & [5] respectively). In these galactic compo- 
nents, one can encounter stars with slightly different [a/Fe]- 
enrichments. For that purpose, we have computed spectra with 
similar atmospheric parameters (T e ff, log(g), [M/H]) and differ- 
ent a-enrichments. We point out that in all these comparisons the 
considered model atmosphere has an [a/Fe] ratio that is consis- 
tent with the one adopted in the synthetic spectrum calculations 
(this is not the case for every AMBRE spectra, see Sect. [5] for 
a discussion of this assumption). As we have mentioned already 
in Sect 2.2, the atmospheric temperature structure of cool dwarfs 



and giants is strongly affected when a-enhancement are consid- 
ered (see Gustafsso n et al. 20081 . 

For instance, Fig. [4] compares the spectra of a metal-rich 
([M/H]= 0.0) giant of the bulge with [a/Fe] =+0.4 and of a typi- 
cal disk giant ([a/Fe]= +0.0) that could be detected in the same 
line-of-sight. Such variations in a-enrichments can also be en- 
countered in bulge giants observed at different galactic latitudes 
(Gonzalez et al., 2011). It can be seen that differences as large 
as 10-20% in relative flux are found for the atomic lines and 
most molecular lines. These differences affect almost any spec- 
tral domain, even those dominated by non-a-rich molecules. 
Furthermore, metallic lines (except those of a-elements) and 
FeH, NH or carbon-rich molecular lines appear stronger in the 
[or/Fe]= +0.0 spectrum. In other words, several lines of non a el- 
ements appear weaker in the a-enhanced spectra, although their 
adopted abundance is kept constant. This is caused by the fact 
that a-elements (and, particularly, Mg, Si and Ca) are important 
sources of electrons. They therefore strongly contribute to the 
total continuous opacity, dominated by the H ion in the stud- 
ied stars (with respect to HI Rayleigh scattering). Thus, the con- 
tinuum flux level of a-enhanced spectra is lower. For instance, 
in the two spectra of Fig. |4j the continuum flux is 25% higher 
in the blue for the [a/Fe] =+0.0 spectrum with respect to the 
[a/Fe] =+0.4 one (this difference decreases towards longer wave- 
lengths and is less than 2% in the near infrared). In summary, in- 
creasing the abundance of a-elements therefore leads to weaker 
lines of non a-atoms and molecules in normalised spectra. 

Similar effects (caused by identical reasons) can be seen 
in Fig|5] where the differences between two halo giant spec- 
tra ([M/H] = -1.5) are shown: one formed in-situ ([o/Fe] = 
+0.4) and one possibly accreted from a satellite dwarf galaxy 
([a/Fe]=+0.0; as found for instance in Fornax or Sculptor, see 
Tolstoy et al., 2009). Again, it can be interestingly noted that the 
main differences between both spectra are found for the atomic 
lines and those of NH, CH, CN and C2 (their lines appear deeper 
in the [a/Fe]= +0.0 spectrum) whereas the OH, SiH and MgH 
lines appear much stronger for [a/Fe] = +0.4. 

As a consequence of all these discussions, it should be 
pointed out again that the integrated spectrum of galactic stel- 
lar components or stellar populations found in external galax- 
ies containing large proportions of cool main-sequence and 
RGB stars is thus affected by their respective enrichment in a- 
elements. This should be carefully considered when interpreting 
such integrated spectra. 

5. Inconsistent [a/Fe] ratios between the model 
atmospheres and the synthetic spectra 

We recall that several spectra of the AMBRE grid have been 
computed from model atmospheres with inconsistent enrich- 
ments in a-elements with respect to the ones adopted in the 
spectral synthesis (a difference, noted Ala/Fe] hereafter, of 
+0.2 dex or +0.4 dex can be encountered^] This shortcoming 
was adopted in order to compute a spectra grid with several 
variations in [a/Fe] (even if no MARCS models were available) 
in order to be able, for instance, to automatically derive the 
[a/Fe] content of the analysed spectra as required by several 
galactic archaeology projects. In 2009, when we recovered the 
MARCS models, too few models with several [a/Fe] -values for 
a given metallicity were available (only some a-poor models 
were existing at that time). The situation is now improving since 

2 The tests performed in the previous section do not suffer from this 
inconsistency. 
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Fig. 3. Normalized spectra computed with only one line specy (all atoms or individual molecules) included in the synthesis. This 
figure illustrates the line blocking by different species for metal-rich stars with solar cr-enrichment ([M/H]= +0.5 & [cr/Fe]= 0.0). 
Examples are shown for a solar type star (T er r= 5750 K & log(g)= 4.5; top left panel), a cool dwarf star (T e ff= 4250 K & log(g)= 
4.5; top right panel), a very cool dwarf star (r er r= 3000 K & log(g)= 4.5; bottom left panel) and a cool giant star (T er r= 4250 K, 
log(g)= 2.0; bottom right panel). To ease the comparison, molecules are grouped with respect to their carbon-rich (C2, CN and CH) 
or or-rich (TiO, ZrO, MgH, CaH, SiH and OH) nature. 



more a-enhanced or cr-poor models have been made available 
since summer 2010 and many more MARCS models should be 
computed in the near future (as they are required by the new 
galactic surveys as GES or Gaia). When such MARCS models 
will be available, the adopted assumption could be easily given 
up when computing a new grid. 



Therefore, the AMBRE grid contains spectra for which the 
continuous and line opacities were computed from (T, r)-laws 
with different [ff/Fe] values. Since most of the a-elements (as 
Mg, Si and Ca) are important electron-donors, their individ- 
ual abundances directly impact the amount of H~ ions (a ma- 
jor opacity source in cool stars) formed in the different atmo- 
spheric layers and, therefore, the computed atmospheric struc- 
ture. Morover, another effect will result from this assumption 
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since the TiO molecule (composed by two a atoms) strongly 
affects the atmospheric structure of the coolest models (T e s<~ 
4 500 K) as already mentionned in the previous section. In order 
to estimate the consequences of this assumption, we have com- 
puted the differences between some AMBRE spectra and oth- 
ers computed without adopting this assumption, i.e. with con- 
sistent [a/Fe] ratios in the model atmosphere and in the spec- 
trum computation (called exact spectra, hereafter). We show in 
Fig. [6] different examples in which the AMBRE spectra have 
been computed with model atmospheres having an [a/Fe] dif- 
fering by Aa=±0.2 dex with respect to the exact spectra. It can 
be seen that the agreement is very good (differences are lower 
than a few percents in normalised flux at any wavelength) for 
r e ff> 4000 K, whatever the value of the gravity is. The largest 
differences are found below A < 5 000 A where most of the lines 
are encountered. However, differences as large as ~10% over 
the whole optical domain appear for much cooler effective tem- 
peratures where many more lines contribute to the spectra. We 
have also checked that these differences, which are probably pro- 
duced by a continuum level not perfectly estimated in the tested 
AMBRE spectra, do not strongly vary when the stellar metal- 
licity increases. Of course, for metal -poor stars the agreement 
becomes much better. 

The effects of larger differences in [a/Fe] -ratios between 
the adopted atmospheric structure and the spectral synthesis are 
shown in Fig.[7]for a rather metal-poor star and a very metal-rich 
one (for which the studied effects are extreme). For rather hot 
stars (T e g> 4500 K), the spectrum differences are still quite small 
whatever the luminosity and the metallicity are. The differences 
are always lower than a few percents for A > 5 000 A and stay 
below ~10% more in the blue. However, larger differences (up to 



10-15%) are encountered for cooler stars with effective temper- 
ature below ~4 000 K. For even cooler stars (M spectral types), 
the quality of the AMBRE spectra with A[a/Fe]=±0.4 dex be- 
tween the model atmospheres and the spectral synthesis is rather 
poor and they should be used very cautiously (or avoided when 
possible). 

Furthermore, in order to better quantify the differences 
caused by a different adopted [a/Fe] ratio between the model at- 
mospheres and the spectral synthesis, we have computed the dis- 
tance (by ax 2 estimation) between an exact spectrum and several 
AMBRE ones with close stellar parameters. Practically, we con- 
sidered AMBRE spectra having T & g, log(g), [M/H], and [a/Fe] 
values between ± 1 000 K, +0.5 dex, ±0.5 dex and ±0.4dex, re- 
spectively, around the parameters of the exact spectrum. We thus 
computed distances between several hundred of AMBRE spec- 
tra and a given exact one. All the stellar parameter combina- 
tions of Fig. [6]&[7] were considered. Such a test allows a check 
as to whether an algorithm could automatically derive the right 
atmospheric parameters of an input spectrum (an exact one in 
this case) using the AMBRE grid. It has been found that the 
right solution is always recovered if the effective temperature 
is larger than ~4500 K, whatever the [a/Fe] enrichment. For 
cooler stars, when A[a/Fe]=±0.2 dex and T e g< 4000 K, the 
right solution is also perfectly recovered for dwarf stars, while 
giants are found ~ 250 K hotter (one grid step in r e ff, all the 
other parameters being well retrieved) when A[a/Fe]= -0.2 dex 
or with a surface gravity slightly larger (one grid step in log(g), 
all the other parameters being well retrieved) when A [a/Fe] = 
+0.2 dex. For the coolest cases (r e ff< ~3 700 K for dwarfs and 
giants with A[a/Fe] = +0.2 dex or when A[a/Fe]= +0.4 dex and 
r e ff< 4250 K), the closest spectra are always found within one 
step grid in each atmospheric parameters. However, depending 
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on the test cases, any of the four atmospheric parameter can be 
erroneously derived although only T e ff or log(g) are affected in 
most cases. In the worst cases (very cool stars and A[a/Fe] = 
+0.4 dex), two of these parameters can be simultaneously badly 
recovered although always within one step grid. We therefore 
conclude that, even in these extreme cases (very cool stars and/or 
large A[o7Fe]) where the accuracy of the parameter determina- 
tion has decreased, the main stellar properties are still enough 
well recovered in order to classify the observed stars. It has also 
been found that the situation improved if only the reddest part of 
the spectra (A >~5 500 A) is considered since a large amount of 
the atomic and molecular lines are found at shorter wavelengths. 

Finally, we would like to point out that similar warnings have 
to be considered for absolute flux spectra. However, most of 
these differences would be partly or completely hidden for ob- 
served spectra with moderate or low Signal-to-Noise ratios (~ 50 
or lower) and/or very metal-poor stars. 



6. Data access 

One of our goals when computing this grid was to offer it to the 
whole astronomical community. It can therefore be easily recov- 
ered upon request to the authors (laverny@oca.eu). Each spec- 
trum (relative and absolute fluxes) can be obtained individually 
in FITS format (3.6Mb each). The whole grid (16 783 spectra) 
has a total size of ~60Gb. 

Also a large part of this grid has been made publicly avail- 
able through the POLLUX database of synthetic stellar spec- 
tra dPalacios etal. 20101) . Every AMBRE spectrum with 5 000 
< ^eff^ 8 000 K can presently be recovered from the POLLUX 
website (http://pollux.graal.univ-montp2.fr/) that contains 7 940 
AMBRE spectra, i.e. about 50% of the whole grid. 



7. Conclusion 

We have computed, and made publicly available, a very large 
grid of FGKM star synthetic spectra covering the whole opti- 
cal domain. It includes almost any possible combination of the 
atmospheric parameters T e ff, log(g), [M/H] and [a/Fe] (includ- 
ing extreme metal-poor and metal-rich stars). These spectra are 
based on (and are mostly consistent with) the MARCS model 
atmospheres and include the most complete atomic and molec- 
ular linelists available so far. We have shown that about 90% 
of the provided spectra, including those spectra with extreme 
[ff/Fe] values that have been computed with [a/Fe] abundances 
inconsistent with those of the model atmospheres, can be safely 
adopted for any use under these criteria: (i) the considered [a/Fe] 
ratios differ at maximum by ±0.2 dex between the model atmo- 
sphere and the spectral synthesis, (ii) high T e g (>~4 250 K) and 
larger A[a/Fe] values, or (iii) the most metal-poor stars for any 
value of A[a/Fe]. 

However, when A[ar/Fe]=±0.4 dex and T eB < -4250 K, the 
spectra should be considered very cautiously (or avoided when 
possible) although they are still similar enough to the exact ones 
(within one step grid in every parameter dimension) to allow 
spectral classification. 

We are well aware that such 'classical' model atmospheres 
and synthetic spectra can be questioned when they are carefully 
compared to real stars, since the effects of the adopted differ- 
ent assumptions (LTE, ID, static atmospheric structure, varied 
abundances of individual elements,...) could be large for some 
specific cases. The accuracy of the adopted linedata together 
with their systematics should also be questionned. For instance, 
for the coolest spectra of this grid (below ~ 2 800 K), we did 
not consider model atmospheres and opacities caused by the 
formation of dust grains. Such an assumption may at least af- 
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feet the red part of the spectra and/or metal-rich models (see 
Chab rier et al., 200"0] l. However, several of these effects could be 
partly corrected by adopting specific calibrations defined with 
rather well known reference stars (see, for instance, the proce- 
dure adopted within the AMBRE project; Worley et al. 2012 1. 
Furthermore, more advanced models (including NLTE effects, 
3D, hydrodynamics,...) are very far from covering the whole 
range of stellar parameters considered in this work, and there 
exists no other alternative than computing such 'classical' grids. 



However the provided grid consists of high resolution spec- 
tra. They can be very easily adapted to any real observed spectra: 
either by slicing them over the covered wavelength ranges, by 
convolution with a well adapted kernel and/or by adopting a spe- 
cific sampling law (see the application to the analysis of FEROS 
spectra in | Worley et al. 2012] ). We have also preferred not to 
consider any stellar rotational velocities and/or macroturbulence 
profiles in order to ensure a very broad use. This grid is at the 
basis of the whole AMBRE project that consists of the analysis 
of most of the ESO high-resolution archive spectra. All of the 
training phases for the parametrisation algorithms have indeed 
been performed with these AMBRE spectra, adapted in turn to 
each of the four selected spectrographs including their different 
instrument set-ups. This grid has also been adopted for several 
tests of the classification algorithms developed for the Gaia/RVS 
spectra and for the preparation of the Gaia-ESO Survey (selec- 
tion of the GIRAFFE setups, estimates of the uncertainties in 
the derived atmospheric parameters and abundances,...). Finally, 
other applications could easily be considered as, for instance, 
the study of the integrated spectrum of the different Galactic and 
extragalactic components, including those with multiple popula- 



tions, together with their property dependencies with respect to 
the stellar [M/H] and/or [ayFe] distributions. 

Finally, the next step in building similar future large spec- 
tra grids would be to overcome the assumption adopted when 
computing synthetic spectra with extreme [a/Fe] values not 
consistent with the available model atmospheres. Considering 
even more complete molecular linelists (including polyatomic 
molecules) would also be of great interest for the coolest spec- 
tra of the present grid. We also point out that an extension to 
hotter T e ff with ATLAS 12 Kurucz stellar atmosphere models is 
planned for the future, adopting the same linelists (extended to- 
wards higher ionisation stages) and consistent solar abundances 
in order to keep a high consistency between different spectra 
grids covering a much larger range of effective temperatures. 
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